Conductivity of suspended graphene at the Dirac point 
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We study transport properties of clean suspended graphene at the Dirac point. In the absence 
of the electron-electron interaction, the main contribution to resistivity comes from interaction 
with flexural (out-of-plane deformation) phonons. We find that the phonon-limited conductivity 
scales with the temperature as T~ v , where rj is the critical exponent (equal to « 0.7 according 
to numerical studies) describing renormalization of the flexural phonon correlation functions due 
to anharmonic coupling with the in-plane phonons. The electron-electron interaction induces an 
additional scattering mechanism and also affects the electron-phonon scattering by screening the 
deformation potential. We demonstrate that the combined effect of both interactions results in 
a conductivity that can be expressed as a dimensionless function of two temperature-dependent 
dimensionless constants, G[T] and G e [T], which characterize the strength of electron-phonon and 
electron-electron interactions, respectively. We also discuss the behavior of conductivity away from 
the Dirac point as well as the role of the impurity potential and compare our predictions with 
available experimental data. 

PACS numbers: 72.80.Vp, 73.23.Ad, 73.63.Bd 
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I. INTRODUCTION 

The discovery of graphene, a single monolayer of 
graphite has initiated a remarkably intensive study 
of electronic properties of graphene structures (for re- 
view, see Refs. |4jl5|). This interest has both fundamen- 
tal reasons and application-related motivations. From 
the fundamental point of view, the interest to graphene 
is largely motivated by the quasirelativistic character 
of its spectrum: charge carriers in graphene are two- 
dimensional (2D) massless Dirac fermions. This leads 
to a variety of remarkable phenomena related to inher- 
ent topology of Dirac fermions as well as to their prop- 
erties in the presence of various types of disorder and 
interactions. Further, the Dirac character of spectrum 
makes graphene a unique example of a system where es- 
sentially quantum phenomena such as the quantum Hall 
effect can be observed up to the room temperature* 6 - 
From the prospective of applications, the technological 
breakthrough in fabrication of flat single-layer 2D sys- 
tems opens a wide avenue for creation of ultimately thin 
2D nanostructures, thus being in the mainstream of the 
general tendency to miniaturization of electronic devices. 
Moreover, suspended graphene samples demonstrate the 
room-temperature mobility as high as 1.2 x 10 5 cm 2 /Vs, 
which is higher than for conventional semiconductor 2D 
structures. Therefore, high-quality suspended graphene 
flakes with the size of the order of 1 /im may show bal- 
listic transport up to the room temperature; 7-15 It is 
widely believed that, in combination with carbon nan- 
otubes, graphene may form a basis for the future carbon 
electronics. Hence investigation of transport properties 
of graphene is a highly topical problem. 

At low temperatures, the resistivity of graphene is 



dominated by scattering off impurities. Away from the 
Dirac point, the dependence of graphene conductivity 
on electron concentration depends strongly on the na- 
ture of scatterersii^ The experimentally observed (ap- 
proximately linear) dependence in most of the samples 
may be explained by strong impurities creating reso- 
nances near the Dirac point ( "midgap states")) 16 ' 17 yield- 
ing a cx nln 2 n, or, alternatively, by Coulomb impurities 
and/or ripples, leading to a cx ni 16 ' 18 " — The dominant 
type (or types) of disorder and the corresponding dis- 
order strength depend, of course, on technology of the 
sample preparation. 

A hallmark of Dirac nature of carriers in graphene is 
the minimal conductivity ~ e 2 /h at the Dirac pointi^ 
Remarkably, it was found experimentally that the mini- 
mal conductivity stays almost unchanged up to a very low 
temperature (~ 30 mK, i.e., three order of magnitudes 
below the impurity-induced transport relaxation rate) i=i 
This can be explained by "protection" of disordered Dirac 
fermions from quantum localization in the absence of in- 
tervalley scattering 2 ^ or in the case of a chiral-symmetric 
disorde n 16 i 23 

At higher temperatures, the graphene resistivity is 
expected to be determined by electron-phonon and 
electron-electron interactions. Manifestations of both in- 
teractions in transport properties o f graph ene have been 
studied in the literature; see Refs. |2J-[34| for discussion 
of the role of electron-electron collisions and Refs. [3oT - l5ll 
for discussion of electron-phonon scattering. In this con- 
nection, two important features distinguishing graphene 
from conventional 2D semiconductor systems should be 
emphasized. First, at the Dirac point of graphene, 
the electron-electron scattering leads to velocity relax- 
ation (though total momentum is conserved just as in 



the conventional case) and, therefore, gives a contribu- 
tion to resistivity.— ~ 27 ' 31 ' 32 Second, a suspended flake of 
graphene is a crystalline membrane, which implies ex- 
istence of specific type of the phonon modes, so-called 
flexural phonons i^ 2 - 

Apart from their role as one of the most important 
scattering mechanisms for electrons, the flexural phonons 
are very interesting from the point of view of mechanical 
properties and thermodynamic stabilit y 53 ' 54 of graphene 
membrane. The out-of-plane fluctuations represent a 
particularly soft mode (oj ~ q 2 dispersion when an- 
harmonicity is neglected versus w ~ q for conventional 
phonon modes), so that they might be expected to be 
very efficient in inducing strong thermal out-of-plane fluc- 
tuations and thus driving the membrane into the so- 
called crumpled phase. This question was intensively dis- 
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cussed in the literature two decades agi 
tion with biological membranes, polymerized layers, and 
some inorganic surfaces (see also the review in Ref. [52| 
as well as more recent papers, Refs. and [HI). It was 
found that anharmonic coupling of in-plane and out-of- 
plane phonons stabilizes the membrane for sufficiently 
low temperatures T, so that the membrane is in the flat 
phase at relatively low T and undergoes the crumpling 
transition with increasing T. The main dimensionless pa- 
rameter characterizing the state of the membrane is the 
ratio of the bending rigidity k to the temperature. For 
graphene, this ratio for room temperature is quite large, 
k/T ~ 30. This reflects a remarkable rigidity of graphene 
and implies that graphene remains in the flat phase up 
to the temperatures several times higher than the room 
temperature. 

In this paper, we explore transport properties of clean 
suspended graphene, with a particular focus on the case 
of zero chemical potential (Dirac point). We show that, 
despite a high bending rigidity, flexural phonons play 
central role in determining the graphene resistivity in a 
broad range of temperatures. We also demonstrate that 
the anharmonicity crucially affects the magnitude and 
the temperature dependence of the resistivity. 

The structure of the paper is as follows. In Sec. [Ill 
we carry out a general analysis of electron-phonon scat- 
tering in a suspended graphene. First, we neglect the 
electron-electron interaction. In this case, the main con- 
tribution to resistivity comes from scattering by defor- 
mation flexural phonons, while other types of electron- 
phonon interaction can be neglected for realistic values 
of temperatures. We predict a power-law dependence of 
the transport scattering rate on the energy and temper- 
ature, l/r tr oc \e\ 2 ' n ~ 1 T 2 ~ ri , where energy, e, is counted 
from the Dirac point and 77 ps 0.7 is the critical expo- 
nent describing renormalization of the flexural-phonon 
correlation functions due to anharmonic coupling with 
the in-plane phonons. As a result, the phonon-limited 
conductivity at the Dirac point scales with the tempera- 
ture as T~ v . 

In Sec. IIII1 we take the electron-electron interaction 
into consideration and demonstrate that its effect is 



twofold: (i) it screens the deformation potential and (ii) it 
induces an additional scattering channel. We show that 
both these effects may be taken into account on equal 
footing by introducing two dimensionless temperature- 
dependent constants, G(T) and G e (T), which character- 
ize the strength of electron-phonon and electron-electron 
interaction, respectively. We find the conductivity and 
demonstrate that it can be expressed as a dimension- 
less function of G and G e (i.e., all the temperature de- 
pendence can be absorbed in these two parameters). 
Different transport regimes correspond thus to different 
parts of (G,G e ) plane (see Fig. |2J). At low temperatures 
electron-electron collisions are more intensive (G e > G), 
while at relatively high temperatures, G becomes larger 
than G e and flexural phonons dominate. When the impu- 
rity scattering is also taken into account, the temperature 
dependence of conductivity at the Dirac point becomes 
strongly non-monotonous at low T. 

In Sec. IIV1 the behavior of conductivity away from the 
Dirac point (non-zero chemical potential /z) is discussed. 
We show that at fi 7^ the interplay of electron-electron 
and electron-phonon interactions leads to a highly non- 
trivial temperature dependence of conductivity contain- 
ing several regions with different power-law behavior. 
Taking into account impurity scattering makes the whole 
picture even more complex. We analyze how the temper- 
ature dependence of conductivity evolves with increasing 
fj, from a non-monotonous (with a maximum in the in- 
termediate temperature range) at low /i to monotonously 
decaying at sufficiently high \x. 

In Sec.[V]we compare our findings with available exper- 
imental data and find a very good qualitative agreement 
between theory and experiment. 

Section PvTl contains a summary of obtained results. We 
also discuss there open questions and prospective direc- 
tions for future research. 

Technical details of the calculation of scattering rates 
and the hydrodynamic approach are relegated to Ap- 
pendixes fA] and [B] respectively. 



II. FLEXURAL PHONONS 

A. Elasticity of graphene 

We start with a discussion of elastic properties of 
graphene. The energy of the in-plane and out-of-plane 
elastic deformations of graphene is given by conventional 
expression known from membrane physics ££L 
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Here u(r) and h(r) are the in-plane and out-of-plane dis- 
tortions, 



Uij = -[diUj + djUi + (dih)(djh)] (2) 
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is the strain tensor, p ~ 7.6 x kg/m 2 is the mass 

density of graphene, A ~ 3 eV/A 2 and /i ~ 3 eV/A 2 are 
in-plane elastic constants, and x ~ 1 eV is the bending 
rigidity. 

As follows from Eq. (JTJ , there are three types of acous- 
tic phonons in graphene: longitudinal (||) and transverse 
(_L) in-plane modes, and out-of-plane flexural mode. The 
mode frequencies read 

W||q = S||<?, OJj_q = S±q, ujq = Dq 2 , (3) 

where s\\ = [(2fi + A) /p] 1/2 ~ 2 x 10 6 cm/s, s x = 

(/i/p) 1/2 ~ 1.3 x 10 6 cm/s, and D = - 0.46 x 

10" 2 cm 2 /s. 

Below, we assume that qa <C 1 (here a is the bond 
length of the honeycomb graphene lattice), thus neglect- 
ing intervalley transitions. In this approximation, we can 
focus on the study of the vicinity of one of two equivalent 
Dirac points, writing the electron Hamiltonian as 

H = hvcrk, (4) 

where u~lx 10 8 cm/s is the Fermi velocity, cr is vec- 
tor consisting of Pauli matrices acting in the sublattice 
space, and wave vector k is measured from the Dirac 
point. Eigcnfunctions of this Hamiltonian and corre- 
sponding eigenenergies are given by 

V'ka = e lkr |Xk), e«(k) = ahvk, (5) 
where a = ±, and 

1 / e -^k/2 \ 

lxi) = V2\ae^ J' (6) 

with ip^ denoting the polar angle of the momentum k. 
The electron-phonon interaction Hamiltonian reads 

H e , ph = J dr*ty eiPh *, ( 7 ) 
where $(r) is the electron field operator, 

* f = S 4aV^a. * = Yl °k«V£a ( 8 ) 

k k 

and V e , p h can be separated into the deformation potential 
(diagonal in sublattice space) and the effective gauge field 
(off-diagonal in the sublattice space) ^ 

V e<P h = V + Va = giuu + g 2 <rA. (9) 

Here g\ ~ 30 eV is the bare (unscreened) deforma- 
tion coupling constant, and g 2 ~ 1.5 eV corresponds 
to coupling to the phonons via the effective gauge field 
A = (A x , A y ) with the components 

A x — QlLxy , Ay — Uxx ^yy • (-^0) 

Equations © and pH|) represent contributions of leading 
order (in gradients of the distortion fields) to the phonon- 
induced scalar and vector potential, respectively. Higher- 
order terms 5 would only give small corrections and are 
not considered below. 



B. Quasielastic scattering by flexural phonons 

Since g\ 3> 52, we will first consider the deformation 
part of the electron-phonon potential that is expected 
to give a dominant contribution to the scattering rate. 
Contribution of the gauge field will be discussed later; 
we will see that it is indeed much smaller. Further, we 
will focus on the contribution of the flexural phonons to 
the deformation potential 

V = 9l (\7h) 2 /2, (11) 

which is larger than that of longitudinal phonons because 
of softer dispersion of flexural modes. 

The transverse displacement field h(r) can be written 

as 




where S is the sample area. Since V is quadratic with 
respect to h, electron in any scattering act emits (ab- 
sorbs) two phonons with the wave vectors qi and q2 , the 
total transferred wave vector being Q = ±qi ± q2 . The 
transport scattering rate, l/r tr (e), can be presented as 
integral over d 2 qid 2 q2- The main contribution to the in- 
tegral comes from the region in (qi, q2) space, where one 
of the momenta, say q2, is much smaller than the other, 
and the integral is logarithmically divergent oc J dqijq?,. 
The upper limit of the integral is given by e/Kv, while the 
lower limit is given by inverse sample size 1/L provided 
that nonlinear interaction of flexural phonons with in- 
plane phonons is neglected. In fact, such anharmonicity 
provides an infrared cutoff for the logarithmic divergency 
due to screening of flexural phonons by in-plane ones^ 7 - 
we will first neglect it and include later into considera- 
tion. Since q 2 -C qi, we find that total transferred mo- 
mentum is given by qi and total energy gained (or lost) 
by electron is given by fio; qi = Dq\ oc hD(e/hv) 2 . This 
energy is much smaller than T [see Eqs. (fTTj|) . (|2TJ)) be- 
low and discussion after Eq. (|20l) ]. so that phonon Planck 
numbers are large and one can replace 6 q in Eq. (1121) with 
^T/hajq exp(— ifq), where y q are random phases (with 
the correlation function (exp[i(ipq — Vq')]) = ^q q')i over 
which the final expression for the scattering rate should 
be averaged. The inequalities fiw qi -C T, ftiOq 2 -C T also 
imply that the phonon potential is quasistatic. Thus, to 
the leading approximation, one can assume that h(r) is 
a static field, 

/ 2T 

H^^^d^TgCOsiqr + ip^). (13) 

q V H ^ 

As seen from Eq. (fT5|l , the rms thermal fluctuation of the 
out-of-plane amplitude, 

^^«\/I7?^ l ' (i4) 
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is proportional to the system size L (when anharmonicity 
is neglected) and the ratio \J (h 2 (r)) / L is controlled by 
the dimensionless parameter T/n. 

Evaluating the transport rate for the scattering on the 
quasistatic random potential given by Eqs. (fTTj) and (fT3|) 
(see Appendix I A we find 



i _ 2 r T2 ln( . 

T tr (e) 7rh\e\ \vL 



(15) 



for a particle with energy e counted from the Dirac point. 
Here 



fJ 



.91 



32x 



5.3 



(16) 



is the dimensionless coupling constant. We see that the 
coupling is quite strong and is additionally enhanced by 
a divergent logarithm. 

As has been already mentioned, the above analysis 
in fact overestimates the thermal fluctuations related to 
flexural phonons. The divergent logarithm in Eq. (|15[) 
appeared because we studied flexural phonons in the har- 
monic approximation. We are now going to take into 
account the anharmonic phonon-phonon interaction. As 
known from the membrane theory ; 52 i anharmonic cou- 
pling of the flexural phonons with the in-plane ones leads 
to screening of the phonon Green function 



G q = (h^) 



(17) 



where hq = bq^/h/2pujq. Specifically, while in the har- 
monic approximation, G q oc T/xq 4 [see Eqs. f| 13[) and 
(|14|)]. the anharmonic coupling suppresses flexural os- 
cillations by modifying the power-law behavior at large 
scales, i.e., for q smaller than certain q c . The value of 
q c for graphene can be estimated by accounting of in- 
teraction between out-plane and in-plane modes in the 
framework of the perturbation theory; 40 ' 51 



9c = 



hv 



= '^-t^ * 18.7 eV. (18) 
4vr^ 2 (2^ + A) v ; 



Let us compare q c to another important momentum scale 
g* that is determined by the condition fiw q _ t = hDq 2 = T, 



9* 



hv 



hv 2 

A, = — ~ 1.25 • 10 3 cV. (19) 



For 91 , 92 <C 9* the electron-phonon scattering is 
quasielastic. Using above estimates for A c and A*, we 
find that 



qr -^qc-^ q* 



(20) 



where qT — T/hv. In the Dirac point, the characteristic 
momentum transferred in a scattering act is of the or- 
der of qx and therefore is small compared to q c . Hence 
the anharmoinic interaction of flexural phonons with in- 
plane phonons should be taken into account. Further, 



Eq. (I20p ensures that relevant momenta are small com- 
pared to , so that the quasielastic approximation used 
above is justified. It is worth noting that this approxi- 
mation also applies away from the Dirac point because 
inequality e/hv < g» (e <C VTA7) is typically satisfied 
for relevant energies and not too small T. 

Interaction between flexural and in-plane phonons 
leads to a power-law renormalization of the bending 
rigidity ! 52 ! 57 ! 5 ^ 



— , for q < q c , 

q 



and the phonon Green function takes a form 



G q — Z 

m \9c 



9 < 9c 



(21) 



(22) 



Here r\ is a critical index and Z ~ 1. On the analyti- 
cal level, the value of r\ was found in the limit of large 
spatial dimensionality^ extrapolation of this result to 
the situation of interest (2D membrane embedded in a 
3D space) yields r] — 2/3. Some modified versions of 
the large-dimensionality approximation have been devel- 
oped, such as the self-consistent screening approxima- 
tion (SCSA) 52 i 57 i 58 and the "non-perturbative renormal- 
ization group" j 5 - 9 - the corresponding results after extrap- 
olation to the physical dimensionality yield r\ = 0.821 
and T) = 0.849, respectively. Clearly, the extrapolation 
is not controlled parametrically; the scattering between 
the above three values may serve as a rough estimate 
of their accuracy. Numerical simulations of the prob- 
lem gave values rj = 0.60 ± 0.10 and r\ = 0.72 ± 0.04 
(see Ref. HH and Ref. HH, respectively). We will use the 
latter value for estimates below ^ As to the numerical 
prefactor Z ~ 1, we did not find its reliable value in 
the literature. In order to obtain Z theoretically, one 
should perform a microscopic modeling of elastic prop- 
erties of a graphene membrane. Alternatively, when the 
anharmonicity regime will be identified in experiment, 
one should be able to find Z from a comparison of exper- 
imental data with the theory. Numerical solution of the 
SCSA equatio n 48 ' 49 yielded Z « 3.5, which is the only 
numerical value available in the literature. It is worth 
noting, however, that SCSA, which becomes exact in the 
limit of infinite dimensionality, is an uncontrolled approx- 
imation for a 2D membrane in a 3D space. For estimates 
below, we use Z ~ 2, which allows us to get a qualitative 
agreement with experiment (see, Sec. W} . 

Physically, the renormalization-induced increase of the 
bending rigidity, Eq. (1211) . is a manifestation of the ten- 
dency of the membrane towards the fiat phase, which is 
realized when the bare rigidity is large: >c/T 3> 1. In 
the opposite situation, the membrane is in the crumpled 
phase; see Ref. [52| for review of the crumpling transition 
between these two phases. As has been mentioned above, 
the ratio x/T is on the order of 30 for graphene at room 
temperature, so that graphene is in the flat phase in the 
whole range of temperatures under interest. 
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Using Eq. ([22)1 , we find that for qi,Q2 <C q c one should 
introduce the cutoff factors Z (qi/q c ) rt and Z (q2/q c ) r> , in 
the integrand in Eq. (|A3j) . Assuming that e <C ^/TA*, 
after some algebra [see Appendix I A 2j we get 



2rp2 



1 2g 2 T 
r tr (e) 7rfi|e| 



(23) 



for lei > VTA; 



CZ 



2 ( jfj ^ 



, for lei « x/TA; 



where C ~ 2.26. We see that for high energies, |e| 3> 
y/TA c , the only effect of the anharmonicity is the re- 
placement of the infrared cutoff, g m i n oc l/L, in Eq. (fT5")) 
with the size- independent value g c . In contrast, at low 
energies, |e| <C \/TA c , the scattering rate is strongly sup- 
pressed and goes to zero as |e| 27?_1 with decreasing the 
energy. The latter case is realized in the Dirac point, 
where |e| ~ T < *frK~ c . 

Equation (|2"3")l will be used below for calculation of the 
Drude conductivity: 



,„„ = e *N I dep(e)(-^) (24) 
de _ .l £|Ttr(e) , , (25) 

*> = sfc (26) 

is the density of states in a single valley, N = 4 is 
the spin- valley degeneracy of the graphene, nj?(e) = 
{1 + exp[(e — [i)/T]} 1 is the Fermi-Dirac distribution 
function, and [i is the chemical potential. The cases fi = 
(Dirac point) and will be discussed in the next two 

sections. 



dni 



v 2 r t r(e) 



e 2 N 



'Tcosli [(e-/i)/2T] 



Here 



III. CONDUCTIVITY AT THE DIRAC POINT 

0* = o) 

A. Interplay of electron-electron and 
electron-phonon scattering rates 

We apply now the above results to the case of zero 
chemical potential (Dirac point). It is instructive to 
begin with substituting the expression for the trans- 
port scattering rate with neglected anharmonic coupling, 
Eq. (|T5j) . into the Drude formula (|2~4"1) . A simple calcula- 
tion yields 



Tt 2 N 



Cph 



h 24g 2 In (q T L) 



(anharmonicity is neglected). 

(27) 



Let us compare Eq. (|27| with the Drude conductivity at 
the Dirac point limited by electron-electron interaction: 24 



Nln 2 



where 



9t 



h 2tt 5 2 (0.69 N — 0.24) : 



l + ( 5 o/4)ln(A/T) 



(28) 



(29) 



is a renormalized constant of electron-electron 
interaction^ g° = e 2 /hnvF is the bare constant, 
A is the ultraviolet cutoff (of the order of the band- 
width), and k is the dielectric constant. In Eq. (|2"9")l it 
is assumed that j° < 1, so that g e remains small in 
course of the renormalization. A more general approach 
presented in Ref. uses 1/N expansion and allows one 
to find renormalization group equations for arbitrary 
<7g, in particular, in the presence of different types of 
disorder (see also Ref. Hi] for analysis of renormalization 
in the presence of optical phonons). 

For graphene suspended in the air, g® is estimated the- 
oretically as g® ~ 2. (In fact, g e may be suppressed by 
using graphene suspended in the media with high dielec- 
tric constant, for example, in conventional water.) In 
the absence of disorder and phonons, renormalization re- 
duces the coupling, so that g e becomes small at suffi- 
ciently low energies (temperatures). Generically, disorder 
and phonons slow down such decrease. Moreover, rela- 
tively strong disorder and/or strong coupling to optical 
phonons may even lead to a non-monotonic dependence 
of g e on Ti 26 ' 41 We do not discuss this case here^ 

Below, we will use the renormalized value of the 
electron-electron coupling, g ei as a parameter of the the- 
ory which can be small, g e -C 1 , or on the order of unity. 
We will also assume the number of "flavors" N to be 
large, N 3> 1. Another important parameter of the the- 
ory, which can be small or large, is g e N (both cases will 
be discussed below). Equation (|28l) was derived under 
the assumption g e N <C 1. In the limit N ^> 1, this equa- 
tion becomes 



In 2 2 



hg 2 0.69 



(30) 



With increasing g e , a oc decreases and saturates (see dis- 
cussion at the end of Sec. lIIIBl) for g e N ~ 1 at the valued 



e 2 A 2 



(31) 



This equation yields an estimate of minimal value of con- 
ductivity in the only presence of electron-electron colli- 
sions. 

Estimating now the conductivity limited by electron- 
phonon and electron-electron interactions according to 
Eqs. ([27)1. (|3"0"j) . and ([21]), we find that the former is 
much smaller, crph/c 00 ~ 10~ 2 . Hence, in the consid- 
ered approximation (neglecting the anharmonicity), the 
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electron-phonon interaction strongly dominates over the 
electron-electron one. 

As discussed in Sec. Ill BL it is important to take 
into account the anharmonic coupling between flexural 
and in-plane phonons that enhances the bending rigid- 
ity and therefore suppresses the electron-phonon scat- 
tering rate. Using the corresponding expression for the 
scattering rate, Eq. ([23| . and taking into account that 
|e| ~ T < \/TA c , we obtain from Eq. ([H} the phonon- 
limited Drude conductivity in the Dirac point, 



Nd 



h 16Z 2 Cg 2 



T 



(32) 



Here, C x = dxx 2 ~ 2tl /cosh 2 (x/2) ~ 2.19. For room 
temperature and r\ = 0.72 (and approximating Z by 2 
as discussed above), we find (A c /T) v ~ 10 2 , so that 
CTph ~ 1.2 e 2 /h. Comparing this estimate with the one 
given by Eq. ([31]) , we see that at room temperature con- 
tribution of the electron-phonon scattering strongly dom- 
inates over the electron-electron one, even for the strong 
interacting case, g e N ^> 1. However, as we will see be- 
low, Eq. (|3"2")l overestimates the contribution of flexural 
phonon scattering because it does not take into account 
screening of the deformation potential. Such a screen- 
ing leads to suppression of deformation potential. As 
a result, contributions of electron-phonon and electron- 
electron becomes of the same order yielding a more re- 
alistic value of conductivity cr oc+p h — 10 ~ 20 e 2 /h (see 
discussion in Sec. IVf . 

Next, we discuss the competition between the electron- 
electron and electron-phonon interactions in the Dirac 
point in more detail. We start with considering a weak 
electron-electron interaction g e N <C 1- This inequality 
ensures that the rate of energy relaxation l/r|f caused 
by electron-electron collisions is much higher than the 
rate of velocity relaxation due to these collisions (in the 
estimates below we assume e ~ T)-M- 



1 



> 



1 



(33) 



Let us first assume that the energy relaxation is also 
faster than the phonon-induced transport rate, 



1 



> 



1 

~ph ' 



(34) 



the relation between 1 / T t e r e and 1 / r tr h being arbitrary [in 
the opposite limit, l/rjf <C 1/T tr h , the electron-electron 
interaction can be neglected and conductivity is given by 
Eq. (|32]l ]. In this first step in calculation of 

conductivity, one should average the inverse scattering 
rate over energy for a given direction of velocity:^ 2 . 



1/T t ° r c ( e ) + 1/T t p r h (e) 



(35) 



where 



(•••)=-/ de(---)p(e)dn F /de 



On the second stage, one should substitute rj 
Eq. (gU). The result reads 



(36) 



tr vo into 



-»CC + 



ph = e 2 N 



e 2 7VTr t a r vo ln2 
2^h 2 



(37) 



Substituting Eq. ([23]) into Eq. (J35J) and using Eq. (76) of 
Ref. H we find 



1 



— ave 
r tr 



T 

Ti 



C 2 Z 2 g 2 



T 



C 3 g 2 e N 



(38) 



where the numerical factors are Ci = 
(2G / tt)[ / dxx 2r > I cosh 2 (x I '2)\ / '[Jdxx/ cosh 2 (x/ '2)] ~ 2.0 
and C 3 ~ 0.989. Using Eqs. §7$ and ||3SJ(, we finally 
arrive at the following result: 



A^ln2 



(J, 



cc+ph 



h 2n [C 2 Z 2 g 2 (T/A c )>i + C 3 g 2 N] 
Equation (|39[) may be rewritten as 



(T, 



cc+ph = 7— r N 2 In 2 S(G,G e 
2nn 



where 



E(G,G e 



1 



G ~\~ G c 



(39) 



(40) 



(41) 



is a dimensionless function of two dimensionless 
temperature-dependent coupling constants, G = G[T] 
and G e = G e [T], defined as 



G e = C 3 g 2 e N 2 



(42) 



(temperature dependence of G e is determined by renor- 
malization of g e ). Let us emphasize that all temperature 
dependence of the conductivity is absorbed into these 
two coupling constants. At low temperatures G <C G e , 
so that we have E = 1/G e and conductivity is limited 
by the electro-electron collisions, while for high tempera- 
tures G 3> G e , the electron-phonon scattering dominates 
andE = l/G. 



B. Screening 

Up to now we neglected screening of the electron- 
phonon interaction by the electron-electron one. It is 
known, however, that such a screening may reduce de- 
formation coupling constant and thus the phonon con- 
tribution to the resistivit y 37 ' 44 ' 45 Quite analogously to 
Ref. |45| we find that Thomas-Fermi screening modifies 
the electron-phonon coupling: 
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1 + 27re 2 NU{Q)/kQ ' 



(43) 



6 



Here Q = ±qi ± q2 is the total transferred momentum 
(see Appendix |X| and n(<5) = U.(oj, Q)\ u =a is the static 
polarization operator. The additional momentum depen- 
dence of the coupling (|4"3")l yields the transport scattering 
rate 



1 



(44) 



1, 



for |e| » ffe iVT, 



calculated in Appendix IA 31 (hereinafter in this section 
for the sake of brevity we omit numerical coefficients on 
the order of unity). Equation (|4"4")) is valid for g e N <C 1. 
In the opposite case, g e N 3> 1, scattering rate is given by 
the upper line of Eq. (|44"]) for all relevant energies, e < T. 

Using dimensionless coupling constants, Eq. (|32]l . di- 
mensional rate 



7 



hN 
Trtt ' 



(45) 



and dimensionless energy x — e/T, one can rewrite 
Eq. flU]) as follows: 



7 ph - G\x\ 



2?)-l 




for |x| < y/G~ e 
for Id > JGl 



(46) 



Since 7 P h oc |a;| 2?7+1 for \x\ — > 0, the conductivity turns 
to infinity due to divergent contribution of small x : 
fph ~ / dx/\x\ 2r > [see Eq. P5]l]. Hence, surprisingly, 
when the screening is taken into account, the electron- 
phonon scattering by itself is not sufficient to yield a finite 
resistivity. Physically, this happens due to enhancement 
of screening with approaching to the Dirac point, which 
results in shunting of the dc current by low-energy elec- 
trons. The divergency is cured by accounting of electron- 
electron collisions whose rate scales at low energies as 32 
7cc ~ At low energies, such collisions win compe- 

tition with electron-phonon scattering even for the case 
G> G e . Consequently, a new low-energy scale, 



G 



2/(4^7+1) 



«r, 



(47) 



appears in the problem, defined by the condition 7 CC ~ 
7 P h- 

Next, we analyze different transport regimes, starting 
from the case G e < 1. To this end, in Fig.[TJ we compare 
7 P h(a;) plotted schematically for G e <C 1 and four dif- 
ferent values of coupling constant G with dimensionless 
electron-electron transport and energy relaxation rates 
[7ee(^) and je(x), respectively] calculated in Ref. [H. 
The averaging procedure appropriate for evaluation of 
the conductivity depends on the relation between je and 
7ph + 7cc- Specifically, for 7^ 3> 7 P h + 7cc one should first 
average the total rate, 7 P h + 7 ee , over energy and use the 



thus obtained averaged rate for calculation of conductiv- 
ity [see Eqs. ([3"5 ]) . (|57 ) . ([55 ]) . and ([55])]. On the contrary, 
for je *C "f P h + 7ee, the conductivity is controlled by 
energy-averaged [see Eq. effective transport time, 

(7ph + 7cc)~ 1 - (In fact, the averaging procedure is only 
important for the numerical coefficient.) 



Ill 




FIG. 1: Schematic plot of the dimensionless electron-phonon 
scattering rate y p h(x) at G e < 1 for four different values of 
G (G increases from I to IV). Also shown are the transport 
scattering rate 7ee(a;) (thick solid) and the energy relaxation 
rate Je(x) (dashed) induced by electron-electron scattering. 

Let us consider regimes I-IV (see Fig. 1) realized with 
increasing electron-phonon effective coupling G: 

• I. G « G e . Electron-phonon coupling is weak and 
7ph 7cc within relevant energy interval (e < 
T, x < 1), so that the phonon contribution to the 
transport rate is negligibly small and S = 1/G e . 
Relevant energies are of the order of temperature, 
x~ 1. 

• II. G e <C G <C Gl~ v . Electron-phonon contribu- 
tion dominates, the screening of the phonons yields 
negligible effect, implying that S = 1/G and x ~ 1. 



III. Gl~ v « G « G 
II, £ = 1/G, x~ 1. 



1/2-2tj 



The same as for regime 



• IV. G X J 2 2r> <C G. The conductivity is determined 
by a competition between electron-electron col- 
lisions and screened electron-phonon interaction. 
The dominant contributions comes from low ener- 
gies, x - (G e /G) 2 / 1+4,) < 1, where 7ph 7oe . 
The dimensionless conductivity S scales with the 
coupling constants as S - (G e /G) 3/1+4 ''. 

Consider now the opposite case G e 3> 1. Calculations 
analogous to the ones carried out in Ref. [12 show that in 
this case, 7 oe ~ \/~x for all relevant energies (x < 1). In 
the absence of phonons (G = 0), the conductivity limited 
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by electron-electron collisions is given by2£ a ~ e 2 N 2 /H, 
and, consequently, E ~ 1. For G^O one should also take 
into account phonons which are strongly screened in this 
case for all relevant energies, so that 7 P h ~ (G/G e )x 2v+1 
for < x < 1. The phonon scattering becomes im- 
portant when G becomes larger than G e . The main 
contribution to the resistivity comes then from the re- 
gion x ~ (G e /G) 2/1+4 '' < 1, where 7 ee ~ 7 P h, yielding 
£~ (G e /G) 3/1+4 ". 



C. Results 

The above results are summarized in Fig. [5J which 
illustrates different scattering regimes in the plane of 
parameters G and G e . In the regions (a) and (b) the 
phonon scattering is weak and the conductivity is limited 
by electron-electron collisions only. Contrary to this, in 
region (c) the electron-electron interaction is weak, the 
phonons dominate transport properties and their screen- 
ing can be neglected. In the region (d), the conductivity 
is determined by competition between electron-electron 
collisions and scattering by screened phonons. As a re- 
sult of this competition, a new energy scale appears in the 
problem, where contributions of both types of scattering 
are of the same order. Finally, on the boundary of the 
region (e) , the conductivity achieves its "quantum limit" 
of the order of e 2 N/h. We expect that the conductivity 
saturates at this value in the whole region (e). 



G e 



(a) 



G e = G 



£ rsj 1 



G "-/ G e = G-*-^ 



x ~ 



/(c) 



J_ 
G 




FIG. 2: Conductivity of clean suspended graphene at the 
Dirac point. Transport regimes characterized by different 
behavior of dimensionless conductivity E are shown in the 
parameter plane of effective (temperature-dependent) dimen- 
sionless couplings G and G e . The analytical expressions for 
dashed lines separating different regimes are given in the white 
boxes. 

It is worth reminding the reader that all the temper- 
ature dependence has been absorbed into dimensionless 
constants G and G e . Since G e depends on T in a very 
slow (logarithmic) manner, the dependence of E (and, 



consequently, a) on T is mostly determined by power- 
law temperature dependence of G. The dependence of 
E(G, G e ) on G for fixed G e is illustrated in Fig. for the 
case of relatively small G e such that (l/N) 21 - 2 ^- 1 )/ 3 < 
G e < 1. These inequalities correspond to a horizontal 
line in (G, G e )-plane (see Fig. [5]) lying above the upper 
left corner of the region (e) but below G e — 1. We see 
that at small G (low temperatures) the electron-electron 
collisions dominate. At intermediate temperatures E is 
limited by phonons (£ ~ 1/G) and at high temperatures 
the conductivity is determined by the narrow region of 
energies where electron-electron collisions and scattering 
on the screened phonons have approximately equal rates. 
In this region, E ~ (G e /G) 3 ^ 1+4 ''. The perturbative cal- 
culations presented above become invalid at very high 
temperature, when E drops down to ~ 1/N and, con- 
sequently, a becomes of the order of the quantum limit 
e 2 N/h. 

In order to make the picture complete, we also showed 
in Fig. [3] the contribution of a static disorder (if ex- 
ists) assuming that it is due to randomly distributed 
charged impurities (another possible type of disorder 
in suspended graphene — adatoms — would lead to simi- 
lar results 1 ^). Such disorder dominates at low temper- 
atures when G is small. Indeed, transport scattering 
rate due to charged impurities is inversely proportional 
to the energy l/r; m (e) oc n.i/e (here m is the impurity 
concentration) and at low T exceeds electron-phonon 
and electron-electron scattering rates at relevant ener- 
gies e ~ T. The impurity-limited conductivity is esti- 
mated as er im ~ (e 2 N/h)TT im (T) cx T 2 /^. This equa- 
tion is valid provided that TT lra {T) 3> H, which implies 
that temperature is not too small. With further lower- 
ing temperature, a lm saturates at the value e 2 N/H (see 
Refs. H6II22II23I for the analysis of the nature of the cor- 
responding quantum critical point and discussion of re- 
lated issues such as localization or antilocalization). As 
illustrated in Fig. [3J impurity scattering becomes rele- 
vant when <7i m becomes smaller than conductivity limited 
by phonons and electron-electron collisions. The decay 
of the conductivity both at low and at high tempera- 
tures can be also understood in the following way. Since 
the phonon potential is quasistatic, it can be treated on 
equal footing with the impurity scattering. Hence one 
could incorporate the impurity scattering into the effec- 
tive coupling constant G, which would then become a 
non-monotonic function of temperature G cx T' n + 1/T 
(we omit temperature- independent coefficients). With 
decreasing temperature, G defined in this way would first 
fall, then reach the minimum, and then start to grow 
again, so that both at very high and very low tempera- 
tures we would arrive at the region (e) ( "quantum limit" ) 
in Fig. dJ 

In the above analysis, we neglected contribution of 
other types of phonons. This can be contrasted with 
the previous publications ; 3 ^ 38 ' 39 ! 43 ' 44 where it was ar- 
gued that the contribution of phonon-induced random 
vector potential might dominate over deformation po- 
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quantum 
limit 




G e N 3 Q 



FIG. 3: Dimensionless conductivity £(G, G e ) at the Dirac 
point as a function of effective phonon coupling G for fixed 
small electron-electron coupling G e . According to the defini- 
tion of G, Eq. (|42p . the shown dependence can be also un- 
derstood as temperature dependence of the resistivity. In the 
low-G (low-temperature) region, the conductivity is governed 
by impurity scattering (shown by dashed line). 



tcntial. Let us estimate the contribution of the gauge 
phonon fields. From Eqs. @ and (ITU)) , one finds the con- 
tribution of the flcxural phonons to the gauge potential: 



A, 



d x hd y h. 



A), 



[{d X hf - (dyhf]/2. 



(48) 



Similar to Eq. (jlip . gauge potential Va is quadratic 
with respect to out-of-plane displacement of the graphene 
membrane. Also, analogously to deformation potential, 
the gauge field is quasistatic. Hence the scattering rate 
may be calculated by using the golden rule for scattering 
on the static potential Va- Proceeding in this way, we 
obtain 



7^ oc G A x 2r >~ 



(49) 



Equation (|49l) differs from the bottom line of Eq. (|4l))) 
only by replacement of G with a much smaller constant: 



G 



2.5 x lQ- d G. 



(50) 



so that the scattering off the gauge field is much less ef- 
ficient than that off the deformation potential. In fact, 
one should be slightly more careful at this point, since, in 
contrast to the deformation potential, the gauge field is 
not screened and Eq. (|4^|) remains valid also at x < \/G e , 
where deformation potential scales as l/G e because of 
screening. Hence, with increasing G e , gauge field may 
come into competition with deformation field. Specifi- 
cally, for 



G e > ( 9 — ) « 4 x 10 2 



(51) 



one could neglect the latter contribution, and the con- 
ductivity would be limited by a combined effect of the 



gauge field and electron-electron interaction. Simple es- 
timates show, however, that the inequality Eq. (|51[) is 
hard to satisfy in realistic system, especially when the 
logarithmic rcnormalization of g e is taken into account. 
Furthermore, there is a second condition: 



G> 



(52) 



which ensures that the rate of electron-electron collisions 
is smaller than the gauge-field scattering rate. Only pro- 
vided that both Eqs. ([51]) and (|52"|) are satisfied (i.e. both 
G and G e are very large), the resistivity is controlled by 
the scattering off the gauge field, E = 1/Ga- This situ- 
ation appears highly unrealistic. We thus conclude that 
the gauge phonon field does not essentially affect trans- 
port properties of graphene in the Dirac point up to very 
high and unrealistic values of G and G e . 

One can also check that in-plane phonons do not give 
essential contribution for realistic values of parameters. 
This is consistent with the previous study^ that came 
to the same conclusion for the case of large chemical po- 
tential fi in the absence of externally applied strain. 



D. Phonon-induced velocity renormalization 

In the previous sections we discussed the electron scat- 
tering rate caused by flexural phonons. Technically, this 
implied calculation of imaginary part of the electron self- 
energy in the quasistatic phonon potential. One may also 
calculate the real part of the self-energy, thus extracting 
information about a phonon-induced spectrum modifica- 
tion. This will allow us to verify the assumption that the 
electron spectrum is not changed essentially, which was 
implicit in our perturbative analysis. 



5v/v 
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nonperturbative 
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f- — i ^ 




Gx 2 " 
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Nx 2(i-n) 



X 



FIG. 4: Velocity renormalization for G e <C 1 and 

G/NG\- V > 1 



Simple calculations yield the following estimate for the 
energy-dependent velocity renormalization, Sv, caused by 
flexural phonons damped by the Thomas-Fermi screen- 
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ing: 



G 



Sv 
v 



r 2 '? 



NG e 
G 1 




(53) 



Comparing Eq. ([53]) with Eq. (f46|) . we see that 

Sv h 
v er tr ' 



(54) 



so that the relative correction to the velocity is on the 
order of the scattering-induced spectrum smearing. In 




FIG. 5: Regions A,B,C where a "non-perturbative" interval 
of energies exists 

most of the cases, spectrum correction is small, Sv/v <C 1 , 
in the whole relevant energy interval x < 1, and thus 
harmless. However, in certain domains of parameters, 
the estimate (f53"|) ceases to be small. This may indicate 
that the calculation of the conductivity within the lowest 
order of the perturbation theory (Born approximation) 
may become insufficient. 

Specifically, for G e > 1 we have Sv/v ~ {G/NG e )x 2 ^ 
at x < 1, so that the maximal value of Sv/v is achieved 
at x ~ 1: <Si; max /u ~ G/NG e - We thus conclude that 
"nonperturbativc" effects (i.e. those going beyond the 
Born approximation) might show up for G e < G/N. 

Consider now the opposite case G e < 1, when the rel- 
ative correction reaches the maximum 



Sv n 



G 



NGl~ v 



(55) 



at x ~ \[G~ e . Hence, "nonperturbative" effects come now 
into play for G e < (G/N) 1 /<- 1 ' r '\ as illustrated in Fig.|U 
From this estimates, we find the regions on the plane 
(G, G e ) for which Sv/v > 1 at some energy interval (but 
at the same time a > e 2 N/K). These regions are marked 
as A,B, and C on Fig. [5] In the regions A and B, the ve- 
locity correction is not small in the interval x\ < x < x 2 , 



where Xt = (G e N/Gf/ 211 and x 2 = {G/N) 1 ^ 1 "^. One 
may expect that this "nonperturbative energy strip" does 
not affect conductivity in the region A, because the width 
of this strip x-i — x\ is much smaller than the tempera- 
ture window: X2 — x\ <C 1. Analogously, the perturbativc 
analysis is expected to give the correct result in region 
C, because Sv/v < 1 in the range x ~ (G e /G) 2/(1+47, \ 
which is determined by the condition 7 p h ~ 7 CC and gov- 
erns the resistivity within the perturbative calculation 
(see previous section). On the other hand, a controlled 
calculation of the conductivity in the region B requires 
going beyond the Born approximation. 

The discussion of such "nonperturbative" phenomena 
is out of scope of the current work and we restrict our- 
selves to a short comment of a somewhat speculative 
character. One may expect that the system becomes 
strongly inhomogeneous, i.e., it can be characterized by 
a local chemical potential /i(r) showing large fluctuations 
around zero (with large number of electrons within each 
such "puddle"). Such fluctuations can be described in 
the framework of random resistor network, which involves 
percolation physics^ Also, within a spatial scale where 
/i(r) is homogeneous, the spectrum of electrons and holes 
might be essentially different from the linear one^i This 
and related issues will be discussed elsewhere. 



IV. AWAY FROM THE DIRAC POINT: 
FINITE n 

In Sec. IIII1 we considered conductivity at the Dirac 
point (p, = 0), which is in the main focus of this paper. 
In the present section, we briefly discuss the behavior 
of conductivity away from the Dirac point, fi =/= 0. The 
phonon-limited resistivity in this regime has been previ- 
ously analyzed in Refs. and H3 where the renormal- 
ization of x was neglected. While this is justified in the 
presence of sufficiently strong externally induced tension, 
the renormalization essentially affects the scattering rate 
when tension is absent (or weak), see Sec. [H] Below we 
explore the effect of flexural phonons on resistivity of 
graphene at nonzero /i in the absence of tension, with 
taking into account the anharmonic renormalization. 

First, we ignore the electron-electron interaction. Sub- 
stituting Eq. into Eq. (gS]), we get 



e 2 N 



(56) 



x < 



16Z 2 C V T 

1 /A c \ v f n \ 2-2-j 



for |^| < T, 



8CZ 2 \T J \T 
1 //i\ 2 I 



for T <H« VTA C 



I 8 \TJ ln(/VVT2£) 



, for VTA C < \n\ 



We see that the conductivity increases with /i as fi 2 2v ~ 
^0.56 fo r T < < y/TA'c and as / In \i for ^JTK C < fi. 
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Now we include the Coulomb interaction into consid- 
eration. Just as in the case fx = (Sec. IHIj) . its role is 
twofold: first, it screens the deformation potential and, 
second, it opens an additional (electron-electron) channel 
of scattering. 

For fi < T [first line in Eqs. ([56])]. the effect of 
the screening was discussed in the previous sections [see 
Eqs. (|M5|) - (|M7]l and For fx > T, the Thomas- 

Fermi screening leads to the renormalization of g (see 
Ref. El) described by Eq. M with 



n(Q) 

Thus, we have to replace 



2tt?i 2 i; 2 ' 



(57) 



(here we omit numerical coefficients of order unity). We 
see that the electron-phonon scattering, being weak at 
fj, = 0, becomes nevertheless dominant at a quite small 
chemical potential fj, ~ T y/ (1 / r ph ) / (1 / t cc ) < T. The 
behavior of the dimensionless conductivity £ with in- 
creasing temperature for /i is shown schematically 
in Fig. [BJ It is assumed in this figure that (i is relatively 
small, so that the temperature-dependent coupling con- 
stant G = G[T] is weak for T ~ /i as compared to the 
coupling constant G e [T] (which depends only logarithmi- 
cally weakly on temperature), 



G[fi] <§C G e [fj] 



(62) 



l + g e Nk F /Q' 



(58) 



where kp = fi/hv. The transport scattering rate is cal- 
culated in Appendix I A 31 For G e <C 1 screening can be 
neglected so that l/r tr and er p h are given by Eqs. (|23|) 
and ([56]) , respectively. For G e 3> 1 , scattering rate reads 



1 2 ff 2 T 2 | £ 



T tr (e) TThGefl 2 



G 3 ln 
CZ 2 



(59) 



( If] x2,? 



for |e|»V7^, 



, for |e|<V7^ 



where C* sa 4. Then, for G e 3> 1, the conductivity limited 
by screened fiexural phonons takes the form 



c p h 



e 2 NG e 




(60) 



, for T<H<CVT^ 



37 G 3 ln(M/VTA^) 



for VtaI< 



Finally, we discuss the role of electron-electron colli- 
sions at ji ^ 0. General equations describing the com- 
petition between electron-electron collisions and electron- 
phonon scattering for arbitrary /i are derived in Appendix 
B. As follows from these results, for // >• T electron- 
electron collisions do not contribute to effective scattering 
rate and the conductivity is given by Eq. (|60p . This im- 
plies that the electron-phonon scattering becomes even 
more important when the chemical potential is tuned 
away from the Dirac point. This conclusion is sup- 
ported by consideration of the conductivity in the region 
< \i <C T for the case, when electron-electron colli- 
sions dominate at fi — 0. As shown in Appendix B, the 
conductivity is given in this case by 



e 2 NT 



1) 



1 



1 



L(l/r ee ) T 2 (l/r ph ) 



(61) 
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FIG. 6: Temperature dependence of dimensionless con- 
ductivity £ for G e < 1 and relatively small /j, such that 
G[/j] <C G e . Dashed lines show conductivity limited by 
charged impurities for different impurity concentration in- 
creasing from the curve (1) to the curve (4). 

For G[T] < G 2 (here G 2 = G[T 2 ] and T 2 » Li should 
be found from the equation G[T 2 ] — G e ^( 2 /T|) the con- 
ductivity is limited by screened phonons. For very small 
G such that G -C Gi = G[n 2 /A c ], the phonons can be 
treated in the harmonic approximation. For G G 2 
the temperature dependence is the same as in Fig. 3. 
At lowest temperatures the conductivity is limited by 
the scattering off static disorder (charged impurities) as 
illustrated in Fig. |S] by dashed lines corresponding to 
different impurity concentrations. Indeed, the phonon 
transport rate, Eq. (|59p . taken at e = [a, has a maxi- 
mum as a function of fj, at /i ~ y/TA c and the maximal 
value decreases with decreasing the temperature as T 3 / 2 . 
Therefore, at sufficiently low temperatures scattering 
by charged impurities with the temperature-independent 
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rate, l/r im oc ni/\e\ 
vant energy interval: 

1 



Ui/fJ,, dominates within the rele- 



> 



(e) r tr (e,T)' 



for \e- fi\< T. 



(63) 



Here Tt r (e, T) is the phonon-induced transport time given 
by Eq. (|2"3")) or Eq. (|59l) ; it is convenient for our purposes 
to write it as a function of two variables e and T. The 
impurity scattering time r; m is in fact also a function of 
two independent variables e and T because of screening 
of the impurity potential. This effect does not change 
qualitatively our results, so that we do not discuss it 
here. 

Let us find the temperature behavior of the impurity- 
dominated conductivity assuming that Eq. (1631) is satis- 
fied. As discussed above, electron-electron collisions be- 
come irrelevant while going away from the Dirac point, 
so that we neglect them and write the conductivity as 

e 2 N , s 

(64) 



a = 



167rfr 2 
dc 



Ttr(e, T)r im (e) 



Tcosh 2 [(e - /i)/2T] r im (e) + r tr (e,T) ' 



Using Eq. ([63]) we expand Eq. (|64|) over power of ri m /r tr 
and keep terms of the zero and first order: 



e 2 N 
16nh 2 



de- 



T im (e) - T 2 m (e)/T tT (e,T)] 
Tcosh 2 [(e-/i)/2T] 



(65) 



Since T <C /i, the integrand is peaked near the region 
e « [i. While calculating contribution of the first term 
in the square brackets we write |e|ri m (e) « £tTi m (jU) + 
[M T im(M)]'( e — A*) + [/ iT im(A t )]"( e — ^) 2 /2- The second term 
in the square brackets of Eq. (|65|) is small and while in- 
tegrating it we neglect temperature broadening of the 
Fermi-function. Doing so, we find for the temperature- 
dependent part of the conductivity, 



5a = a{T) 



a(0) 

TT 2 T 2 
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(66) 



lfmm(p)]" - 



This equation is valid for an arbitrary type of impurity 
scattering. For charged impurities we find that 5a is 
given as a sum of two terms of different signs: 



T 2 /i 3 1 

5<j(T) CX x — r 

ni nf r tr (At,r) 



(67) 



(here we omit coefficients which do not depend on fi, T 
and rii). The analysis of Eq. (|67|) shows that "metal- 
lic" behavior of conductivity (i.e. its increase with low- 
ering T) at low rii changes to an "insulating" one with 
increasing n, as illustrated in Fig. [6] by dashed lines. Al- 
ternatively, this crossover in the behavior of conductivity 
at relatively lower temperatures can be observed if one 
changes fi at fixed disorder. One can also see from Fig. [6] 
that at intermediate values of rij (or of /x) the dependence 
of conductivity on T may have two maxima [see curve (2) 
in Fig. |B] . 



V. COMPARISON TO EXPERIMENT 

We now compare our results with available experimen- 
tal data£ (see also Ref. |47l ). To this end, we plot con- 
ductivity as a function of temperature and chemical po- 
tential (or, equivalently, electron concentration) for the 
same values of parameters as in Ref. H. We assume that 
electron-electron coupling is renormalized from the value 
5g at large energies (~ A) to g e ~ 1 in the room tem- 
perature interval, so that G e — C%g 2 N 2 ~ 16 3> 1. The 
conductivity is obtained by interpolating of equation for 
total scattering rate (which is the sum of the electron- 
electron, electron-phonon, and impurity scattering rates) 
between asymptotical expressions presented in the previ- 
ous sections and by substituting thus found rate into the 
Drude formula Eq. (|24[) . The main purpose of such an 
interpolation is to estimate characteristic values of con- 
ductivity and analyze qualitatively the dependence of a 
on (i, T and rii. Due to evident reasons, we do not pretend 
to get quantitative agreement with experiment. First of 
all, interpolation procedure yields numerical value of con- 
ductivity up to the coefficient on the order of unity at the 
boundaries separating regions of parameters correspond- 
ing to different transport regimes. Second, even asymp- 
totical expressions contain some unknown numerical co- 
efficients. In particular, the coefficient Z ~ 1 entering the 
phonon correlation function [see Eq. (|22[) ] is not known 
as it was discussed above. Below we use Z = 2, which al- 
lows us to get reasonable agreement with experiment. We 
also do not know the numerical coefficient in the equa- 
tion 7 ee ~ sfx for dimensionless electron-electron scatter- 
ing rate in the limit G e ^> 1. In the estimates we choose 
this coefficient to be 1/4 which yields a best fit to ex- 
perimental data. Another issue, which was not resolved 
rigourously, is the screening of the impurity potential. 
In the estimates we simply assume 1/T; m ~ e rii/\e\K 
and choose the coefficient in this equation to be unity. 
This yields a good approximation for impurity scattering 
rate (up to a numerical coefficient) at least for not too 
low temperatures, when conductivity is much larger than 
e 2 /h. 

At the Dirac point the data shown in Fig. 4 of Ref. 
show an increase of conductivity by factor of 2 3 when 
temperature increases from 5 to 200 K, with a clear sat- 
uration around 200 K. Let us compare this picture with 
theoretical estimates. The calculated conductivity at the 
Dirac point is plotted in Fig. [7] for different values of 
impurity concentration rii. A curve corresponding to 
rii = 0.5 x 10 10 cm -2 (third curve counted from the bot- 
tom) reasonably agrees with experiment. A similar be- 
havior with a somewhat larger ratio er(200 K)/cr(5K) is 
seen in Fig. 3 of Ref. H3- We expect that with further 
increase of temperature the conductivity will drop due 
to electron-phonon scattering, as found in Sec. IHII of the 
present work and is shown schematically in our Fig. [3] 
As seen from Fig. the drop of the conductivity with 
T can be observed at smaller T provided that one uses 
cleaner samples. 
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r(K) 

FIG. 7: Conductivity at the Dirac point (fi = 0) for differ- 
ent values of the impurity concentration (n^/lO 10 cm -2 
If) 10" 3 , 10" 2 , 10" 1 , 0.5, 3, 10) increasing from the 
top to the bottom. Dashed line corresponds to SCBA limit 
a — 4e 2 /irh. Calculations are controlled only well above this 
line. 



Let us note that the conductivity curves in Fig. [7] go 
to zero in the limit T — > 0, in view of the vanishing 
density of states at the Dirac point. If we would in- 
clude disorder-induced level broadening self-consistently 
[in the framework of the self-consistent Born approxima- 
tion (SCBA)], we would get instead a limiting conduc- 
tivity value 4e 2 /irh (marked by a horizontal dashed line 
in the plots). The actual behavior of conductivity in this 
regime is controlled by quantum interference effects that 
lead to localization, antilocalization, or quantum criti- 
cality, depending on the character of disorder. 22 In the 
present paper we do not discuss these phenomena, as our 
focus is on regimes where the dimensionless conductivity 
is sufficiently large and quantum interference corrections 
do not change it significantly. 

When one moves away from the Dirac point, the ex- 
perimentally observed temperature dependence changes 
(see Fig. 2 of Ref. HI), and at sufficiently large /j, the con- 
ductivity becomes monotonously decreasing function of 
T. This evolution is in a very good qualitative agree- 
ment with our results. To see this, we fixed the impurity 
concentration at the level = 0.5 x 10 10 cm -2 , and cal- 
culated resistivity (just as in Fig. 2 of Ref. @) as a function 
of the electron concentration n = N/j, 2 /AttH 2 v 2 for differ- 
ent values of temperature (the same as used in Ref. 0) . 
The results are plotted in Fig. [8] and look very similar to 
the ones presented in Fig. 2 of Ref. [H. Within the grey 
area temperature dependence is "insulating" , while out- 
side this region it is "metallic" . One of the main features 
of this picture is the existence of the "stationary" point, 
where T-dependence changes. To illustrate the existence 
of this point in a more transparent way we plotted in 
Figs . 151 and [TUl the conductivity as a function of temper- 
ature for fixed impurity concentration and different val- 
ues of the chemical potential. The "stationary" point in 
Fig. [5] corresponds to existence of more or less horizontal 
lines in Figs. [5J and [TU] separating regions with "metallic" 
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FIG. 8: Resistivity as a function of electron concentration 
at m = 5 x 10 9 cm -2 and different temperatures (T/1K = 
5, 40, 90, 150, 230) increasing from the bottom to the top 
at large n. Within the grey area temperature dependence is 
"insulating" , while outside this region it is "metallic" . 




100 200 300 



T(K) 

FIG. 9: Conductivity at fixed impurity concentration 
(rii — 5 x 10 9 cm -2 ) for different values of chemical poten- 
tial O/IK = 0, 50, 70, 90, 110, 130, 150, 170, 190, 210, 230, 
250, 270, 290, 310, 330, 350, 400, 450, 500, 550, 600) increas- 
ing from the bottom to the top. Dashed line corresponds to 
SCBA limit a = 4e 2 /irh. 




100 200 300 



7(K) 

FIG. 10: Conductivity at fixed impurity concentration (m = 
10 9 cm -2 ) for different values of chemical potential fju/1 K = 
0, 10, 20, 30, 35, 40, 45, 50, 60, 70, 80, 90) increasing from 
the bottom to the top. Dashed line corresponds to SCBA 
limit a — 4e 2 /irh. 



and "insulating" behavior. As seen from Figs. [51 and [TU1 
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the transition between different types of T-dependence 
becomes more pronounced with decreasing the impurity 
concentration. It is worth noting that the transition may 
be also obtained for fixed p by changing the impurity 
concentration (for example, by annealing the sample) as 
illustrated in Fig. [TT] 




l , , , , , 

100 200 300 



r(K) 

FIG. 11: Conductivity at p = 50 K for different 
values of the impurity concentration (n^/lO 10 cm -2 = 
0.07, 0.085, 0.1, 0.125, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.7, 0.9) 
increasing from the top to the bottom. Dashed line corre- 
sponds to SCBA limit a = 4e 2 /nh. 

Sufficiently far from the Dirac point, such that /i>T 
for all relevant temperatures, only the first two phonon- 
controlled regimes of Fig. [5] survive, implying a crossover 
from the conductivity scaling a cx p 2 /T 2 ln(p 2 /T A c ) at 
lower temperatures (T < p 2 /A c ) to a cx p 2 ~ 2v /T 2 ~ v at 
higher temperatures T ^> p 2 / A c . However, the tempera- 
ture, separating two regimes, p 2 / A c , turns out to be very 
small, on the order of 1 K, even for largest chemical po- 
tentials, p ~ 600 K, studied in Ref. @- In other words, ex- 
perimental situation corresponds to anharmonic regime. 
While the data shown in Fig. 3 of Ref. 91 do indicate 
a power-law increase of resistivity with temperature (at 
relatively large temperature) , the exponents do not quite 
agree: the theoretical dependence, p ~ T 2 ~'' ~ T 1,3 [see 
upper line of Eq. (|60[) ] , turns out to be slightly stronger 
than experimentally observed linear one, p ~ T. 

Consider now the dependence of the conductivity on 
the electron density. As observed in Ref. [|| this depen- 
dence is qualitatively different in clean and dirty sam- 
ples. When impurity concentration is large, conductivity 
is a linear function of the concentration, a ~ n, while at 
the same sample after annealing the dependence becomes 
sublinear: a increases with electron density n cx p 2 with 
an exponent considerably smaller than unity, which is at 
least in qualitative agreement with the a cx p 2 ~ 2v cx n 1 "'' 
predicted above for transport away from the Dirac point 
[see upper line of Eq. (foT)|)]. This is illustrated in Fig.lHfl 
where a is plotted as function of n for fixed tempera- 
ture, T — 40 K (the same as in Ref. 0), and two different 
values of the impurity concentration. This picture looks 
very similar to the Fig. 1 in Ref. 0. Physically, the linear 
dependence at large rij is caused by impurity scattering, 



while sublinear one at low is due to the phonon scat- 
tering. 





200 




« ; = 5xl0 9 cmJ/ 




150 














^6 


100 
50 


7"=40K \ 


\ / -1 1A 10 2 

\ / h,= 3x10 cm^ 



-2.5 -1.5 -0.5 0.5 1.5 2.5 

n (10 cm ) 

FIG. 12: Conductivity as a function of electron concentration 
at T — 40 K and two different values of impurity concentra- 
tion: m = 5 x 10 9 cm -2 and m = 3 X 10 10 cm -2 . 

Finally, we make a rough estimates of the numerical 
values of the conductivity. For this purpose, we con- 
sider temperature dependence of the conductivity for 
m = 0.5 x 10 10 cm~ 2 and different values of p (see 
Fig. For p = 600 K (which corresponds to the high- 
est density n e — 2 x 10 11 cm" 2 studied in Ref. [^) and 
T = 200 K (which is approximately the highest tempera- 
ture in Ref. 0) theoretical prediction is a ~ 50e 2 /h (up- 
per curve in Fig. |9|). In the Dirac point theory predicts 
a ~ 9e 2 /h (lowest curve in Fig. The experimental 
data of Ref. i yield a ~ 70e 2 /h for p = 600 K and 
a ~ 12e 2 /h for p = 0, in reasonable agreement with 
our findings. 

It was suggested in Refs. I45ll47l that the experimental 
data of Refs. I9II47I might be to some extent affected by 
strain that may result from fixing the sample at the con- 
tacts. It is worth emphasizing that in the present paper 
we have studied the case of a strain-free graphene and 
obtained reasonable qualitative agreement with experi- 
ment. Furthermore, the strength of the tension depends 
on the procedure of preparation of a sample. In partic- 
ular, the built-in tensions in freely suspended graphene 
monolayers produced by means of chemical reduction of 
graphene oxide were founcU^ to be considerably weaker 
than in mechanically exfoliated graphene samples. 



VI. SUMMARY 

To conclude, we have studied transport in suspended 
clean graphene in a broad range of temperatures. We 
have explored the interplay of electron-phonon and 
electron-electron interactions and have found that the 
scattering off flexural phonons controls the resistivity at 
relatively high T. Taking into account the anharmonic 
coupling of flexural and in-plane phonons was crucial for 
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correct evaluation of the graphene conductivity. 

Our results for the conductivity can be expressed in 
terms of two dimensionless coupling constants G and G e , 
characterizing the strength of the electron-phonon and 
electron-electron scattering, respectively. Both constants 
depend on temperature: G shows a power-law scaling, 
G ~ T v , with the exponent r\ describing the scaling of the 
bending rigidity of graphene membrane with the length 
scale due to anharmonicity, while G e slowly (logarithmi- 
cally) changes with T due to renormalization of the Fermi 
velocity caused by electron-electron interaction. 

At the Dirac point, the dimensionless conductivity 
E(G, G e ) of clean suspended graphene depends on T only 
through temperature dependence of G and G e . At low 
temperatures, G < G e and phonon scattering is negligi- 
ble. In the room temperature interval G> G e and trans- 
port is dominated by phonon scattering. At sufficiently 
high temperatures one should take into account screen- 
ing of the deformation potential due to electron-electron 
interaction. The overall "phase diagram" of scatter- 
ing regimes and conductivity scaling of clean graphene 
at the Dirac point is shown in Fig. [2j A characteris- 
tic temperature dependence of conductivity is shown in 
Fig. 03 There, we also included the regime of low tem- 
peratures where the resistivity is controlled by disorder. 
[The lowest temperatures where quantum criticality or 
localization (antilocalization) effects come into play are 
not considered in this paper, see Refs. [I6II22II23I ] As 
seen in this figure, the conductivity first increases with 
T (disorder-dominated regime), then shows a plateau 
(electron-electron scattering), and then drops due to 
electron-phonon scattering. In the high-temperature part 
of the latter regime, the Thomas-Fermi screening of 
phonons becomes important. Remarkably, at still higher 
temperatures the system enters the ultimate quantum 
regime (conductivity of the order of e 2 /h analogous to 
the one at lowest temperatures). In view of the qua- 
sistatic nature of flexural phonons, quantum interference 
phenomena are expected to be relevant in this regime, 
despite rather high temperatures. 

Away from the Dirac point (/i 7^ 0), the role of the 
electron-phonon interaction increases as compared to the 
electron-electron one (since the latter conserves momen- 
tum). In particular, the electron-phonon scattering, even 
if weak at /x = 0, becomes dominant at a quite small \i 
such that (J, -C T [see Eq. ([61]) ]. The temperature de- 
pendence of conductivity away from the Dirac point is 
sketched in Fig. [6] We also discuss the effect of static 
disorder (assuming for definiteness that charged impuri- 
ties are the main source of disorder) as shown by dashed 
lines in Fig. [SJ With increasing /1 (or, else, with lowering 
impurity concentration at fixed nonzero fi), the behav- 
ior of the Drude conductivity in the disorder-controlled 
low-temperature regime changes from "insulating" (a de- 
creases at lowering T) to "metallic." 

Our findings qualitatively agree with experimental 
data of Refs. 1 and S3, see Sec. and Figs. EHU We 
hope that the results of this paper will stimulate further 



experimental investigations of conductivity of suspended 
graphene, including a systematic investigation of temper- 
ature dependence at different chemical potentials for T 
up to (or even higher than) room temperature. 

A number of problems related to this research have 
been left open. First, this includes the possibility of es- 
sentially non-Born and of "ultimate quantum" regimes 
at high temperatures at the Dirac point. Second, the 
detailed analysis of the scattering by charged impurities 
with the account of self-consistent screening is required. 
It is also interesting to study the effect of flexural phonons 
on quantum transport in suspended graphene in trans- 
verse magnetic fields. 
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Appendix A: Calculation of scattering rate 

In this Appendix, we present a derivation of the trans- 
port time due to scattering off flexural phonons. 



1. Harmonic approximation 

We first neglect the anharmonicity. The quasistatic 
random potential representing approximately the dis- 
placement field of flexural phonons has the form: 

V(r) = — > ^-^-sm(qir + ( / 3 qi )sm(q 2 r + ^ q2 ). 

(Al) 

[Eq. (jAlj) is obtained by substitution of Eq. (fT3| into 
Eq. (jlip .j Averaging squared matrix element of transi- 
tion between ■0k,a and ipk\/3 over the phases y q , we find 
that k' — k ± qi ± q 2 , where all four combinations of + 
and — in front of qi and q 2 should be taken into account: 



9lT 2 ^ (qxq,) 2 a , 



E 



'\,:ix^)l 2 ^,k ±qi±q2 . 

The transport scattering rate for an electron in the 
branch a with the energy e = avk is given by 



\x 2 S 2 ^. q\q\ 
qi,q2,± 1 z 



1 



_ K9jT 2 V f d 2 qi d 2 q 2 3 

r£(e) " 4Sx 2 f^J (2tt) 2 ( 27r )2 IUklx k± qi ±q 2 



(qiQ2)%, , 



— 4-4-^ e k± qi±q2 " ^Xl - n k n k±qi±q2 ). 
I1I2 



(A3) 
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In the quasielastic approximation, a = (3 due to delta- 
function in Eq. (|A3j) and r tJ . is the function of |e| = vk. 
Taking into account that the dominant contribution to 
the integral comes from the region where one of the mo- 
menta is much smaller than another one, say q 2 <C qi , and 
using |(x£lXk')l 2 = (l + n k n k /)/2, we reduce Eq. to 
the form 



1 = TrgfT 2 v f dq 2 d 2 qi 
r tr (e) 8fi^ 2 fW 2^q 2 (2^) 2 g 2 



x ^ ±qi ™e£)[l-(n k n k±qi ) 2 ]. (A4) 
Using the identity 

S(k - |k + q|) = 5{q - 2kcoS(f)/\ cos<^| (A5) 
(here ip is the angle between q and k), we get 
1 



T tr (e) Z2ir 2 h 2 x 2 v J q 2 J qi J_ 7r / 2 



6(qi — 2/ecos tp) 9 
— 1 — cos 



COS if 



(A6) 



Finally, taking also into account contribution of the re- 
gion q 2 ^> qi, we obtain 



1 



rt r (e) I6irh\e\>c 2 



9\T" ^ / 9max 



(A7) 



where the energy e is counted from the Dirac point, 
q m - m oc l/L and g max oc |e|/w are infrared and ultravio- 
let cutoffs, respectively. This yields Eq. ([T5l) of the main 
text. 



2. Including anharmonicity 

Now we include the anharmonicity of flexural phonons 
and derive Eq. (|23p of the main text. In order to treat 
both harmonic and anharmonic regions of momenta, we 
introduce in Eq. (| A3|) a factor ((71)8 ((72), where 



e(«) 



1, for q^> q c 



Z(q/q c ) r i, for g < g c 
We also use the identity 



(A8) 



)r] 



//• d 2 Qd 2 v 
d 2 Q<5(Q± qi ±q 2 ) =y --gp- exp[i(Q±qi±q : 

(A9) 

which allows us to replace ±qi ± q2 with Q in 
£ k ± qi ±q 2 ' n k±qi±q 2 and Xk± qi ± q2 , and integrate in 
Eq. (|A3I) first over d 2 qi<i 2 q2. Denoting the result of in- 



tegration as £(r) we find: 

d 2 qi d 2 q 2 (qiq 2 ) 2 




(2^) 2 (27r) 2 qfql 



1 



e( qi )Q(q 2 )e 



(±qi±q 2 )r 



Jo q 



1 

q c r 



for q c r <C 1 



q c r 



2'/ 



(A10) 



for q c r > 1, 



where £* = [/ °° fte"" 1 J (a)] 2 + [/" fte"" 1 J a (x)] 2 ~ 
2.77 (for 77 = 0.72). Next, we substitute Eq. (lATO)) into 
Eq. (|A3I) and integrate over angle of vector r. We get 



1 



g 2 T 2 



Ttr(e) 



where 



(All) 



/ (^^^+Q- e k)[l-(nkn k+Q ) 2 ]A(Q), 
A(Q) = 4n 2 Q 2 J d 2 re- lQr £(r) (A12) 



7r Q / drrJ (Qr)Z(r). 



By using property of Bessel function xJq(x) — 
d[xJ±(x)]/dx we integrate by part and find asymptotics 
of the function A(Q) 



A(Q) 



(A13) 



8tt 



In ( ^ ]. lor (.) <Ir 



CvZ' 2 ( - 



2</ 



dxx 2n Ji(x), for Q < <? c . 



,9c/ Jo 

Using Eq. (|A5|) we find from Eq. ([ATT]) . 
1 ,g 2 T 2 



47T 3 ?l 2 U 

fdQ r /2 , 5(Q-2kcosip) 2 „ . 

x /-^/ d V— ^(l-cos 2 2<^)A(Q) 

J ^ J-n/2 COS ^ 



2y2 /-Tr/2 



27T 3 ft|f 



7r/2 



d(^sin 2 (^)A(2fccos^). (A14) 



Substituting here asymptotics of A(Q) we finally obtain 
Eq. (|2"3")l with the coefficient C given by the following 
equation: 



C 



2.26. 



/2 
tt/2 



cfea; 2r, J\(x) / c?^ sinV cos 2 V 
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3. Role of the screening 



Appendix B: Hydrodynamic approach 



Here, we calculate the phonon-induced transport scat- 
tering rate in the presence of the screening due to e-e 
interaction. We start with the case of fi -C T. Replacing 
in Eq. (|43j) e 2 /Hud with renormalized value g e and using 
expression for II in the Dirac point obtained in Ref . |32| we 
conclude that screening results in the following replace- 
ment: 



9 



l + g e Nf(Q/2q T )- 



(A15) 



where 32 



/(*) = ~ 



oo 



sinh(uz)w(l — v 2 ) + smh(v z)u(u 2 — 1) 



< J duldv- — . _ 

/i Jo [cosh(uz) + cosh(vz)]uvy / (u 2 — 1)(1 — v 2 ) 



\n.2/z, for z < 1 
7r/8, for z > 1. 



(A16) 



The screening leads to the following modification of the 
expression Eq. (|A14|) for the phonon-induced scattering 
rate: 



1 



g 2 T 2 



1-/2 sin 2 (ip)A(2k cos ip) 
^[l + <? e A/(6Cos^/r)] 2 - 



(A17) 

Combining Eqs. (|A13|) . ([AlB) , and (|AT7)) . we obtain 
Eq. (gl]) of the main text. 

Away from the Dirac point, for /i > T, the Thomas- 
Fermi screening is accounted for by the replacement ([55]) , 
which can be rewritten as 



1 + y/G e k F /y/C 3 Q 

Then Eq. (|A14|) becomes 



(A18) 



g 2 T 2 



r P h ( e ) ir^h\e\ J 



-it/2 



dip 



sm 2 {(p)A{2k cos 



1 



/a 



2\fCzk cosip 



2- (A19) 



We see that screening can be neglected for G e <C 1, while 
for G e > 1 we get 



1 



9 2 T 2 C 3 \e\ 

TT 3 HG e fJ, 2 



tt/2 



dtp sin 2 (2ip)A(2k cos </?) . 



(A20) 

From Eqs. (jAl3|t and (|A20|) we find that for G e > 1 
Eq. (j2"3")l is replaced with Eq. (I5i?]) of the main text, where 



C = CC : 



f-n/2 d<P sm2 cos 2?) ip 



f-!v/2 Sm2 V COS 2 * 7 (jf 



Here we present details of the hydrodynamic approach 
used for calculation of the conductivity in Sec. [IV] (see 
also Refs. H^-HIHllll) . We start from kinetic equa- 
tion (we keep H = 1 throughout calculations restoring it 
in the final equations) 



dn a 



eE 



dn a 



(Bl) 



where I oc and I p h are electron-electron and electron- 
phonon collisions integrals and the index a = ± labels 
bands of positive and negative energies. Let us introduce 
two currents: 



J = / d 2 k n a k, 

a J 

j = ^2 J d 2 k n a v. 



(B2) 
(B3) 



In a conventional semiconductor with quadratic spec- 
trum, these currents are proportional to each other. Im- 
portantly, this is not the case for graphene, so that veloc- 
ity may relax even for a momentum conserving scatter- 
ing such as electron-electron scattering. For simplicity, 
we will assume that 1/te due to electron-electron scat- 
tering is much larger than other scattering rates. The 
peculiarity of the kinematics of particles with linear dis- 
persion yields fast equilibration of carriers within a given 
velocity directioni 24 i 25 ' 27 i 28 Therefore, the electron gas in 
graphene is described by the Fermi distribution function 
characterized by local temperature and chemical poten- 
tial both depending on the velocity angle. Using two 
variables ; 68 ' 69 electron energy and the velocity unit vec- 
tor v = v/v, instead of electron momentum and band 
index a, the distribution function takes the form 



n(e,v) 



exp{[e-M(v)]/T(v)} + l- 



(B4) 



We assume that electric field E is small, so that T(v) — 
T oc Ev, /i(v) - // oc Ev. Expanding Eq. (IB4[) up to 
the first order with respect to E, we find the following 
expression for correction to the Fermi distribution func- 
tion: 



5n 



driF eEv 
~de T 



-X- 



Here 



X = Xo + Xie/T, 



(B5) 



(B6) 



and xoj Xi are energy- independent amplitudes (in a non- 
stationary case, these amplitudes depend on time). The 
currents (|B2[) and (|B3[) may be written as 

J=(J t >, j = (je), (B7) 
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where the current densities in the energy space J e and j e 
are expressed in terms of xo and Xi : 



eET 
eEv 2 



Xo + ■jsXi 



(B8) 
(B9) 



and (■ ■ ■ ) is given by Eq. (1551) . Due to the fast energy 
relaxation one may reduce the kinetic equation to the 
simple balance equations for J and j, or, equivalently, 
to the equations for \o an d Xi- To this end, we multi- 
ply Eq. (IB1[) by k and v and integrate over energy and 
velocity angle taking into account that electron-electron 
collisions conserve momentum, while velocity may relax. 
The result reads (see Refs. 13211681 for discussion of the 
properties of I ee ) 



which directly follows from Eqs. (IB3|) . (IB5|) . and (|B6[) . 

For \x = 0, coefficients (e/T) , (e/Tr ph ) including aver- 
aging of odd functions of energies turn to zero, so that 
we find 



Xo 



(1) 



(1/Tee + 1/Vph) ' 



Xi =0, 



(B15) 



and restore Eq. ([571) for conductivity. 

In the opposite limiting case, \x T, all averages en- 
tering Eqs. (jB12j) and (jB13j) are calculated with the use 
of equation (A(e)) A(p)p(/j,) and we find 



8J eE (e) 

dj _ eEv 2 (l) 
dt 2 



TphO)/ ' 

% - e dje/de 



v T ph(e) ; 



(BIO) 
(Bll) 



By using Eqs. (|B8 |) — (|B11|) one may derive equations de- 
scribing relaxation of the amplitudes xo and xi to their 
stationary values. The latter can be found from the fol- 
lowing set of equations: 



t »^£) + »(k£>- ,BI2i 



(1) 



Xo 



1 



1 

Tph 



Xi 



€ 1 

Tr ph 



(B13) 



The solution of Eqs. (|Bl2|) and (lB13|) should be substi- 
tuted into the expression for conductivity, 



e 2 v 2 N 



(B14) 



T 

Xo = ,xi = T P h(p) — ■ 
A* 



(B16) 



The conductivity is given by 



(j = £!^giM (B17) 

and does not depend on the rate of electron-electron col- 
lisions. Using Eqs. (|59|) and (|B17|) we arrive to Eq. (|60|) 
for the conductivity. 

Consider now the case < /i -C T. Remarkably, a 
new regime arises in this region provided that electron- 
electron collisions dominate over phonon scattering at 
/i = 0. To see this, we first notice that (e/T) ~ 
H/T, (e/Tr ph ) ~ (/i/T)(l/r ph ) and {e 2 /T 2 r ph ) ~ 
(l/r p h), where (l/r p h) is calculated for /i = 0. Then 
from Eqs. (|BT2| - (|Blt)) we find Eq. (f5T|). 
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